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We study a particular nonlinear dispersion relation uj p {k p ) — a series expansion in the physical 
wavenumber k v — for modeling first-order corrections in the equation of motion of a test scalar 
field in a de Sitter spacetime from trans-Planckian physics in cosmology. Using both a numerical 
approach and a semianalytical one, we show that the WKB approximation previously adopted in the 
literature should be used with caution, since it holds only when the comoving wavenumber k S> aH. 
We determine the amplitude and behavior of the corrections on the power spectrum for this test 
field. Furthermore, we consider also a more realistic model of inflation, the power-law model, using 
only a numerical approach to determine the corrections on the power spectrum. 



I. INTRODUCTION 

Inflationary models provide answers to many problems 
in standard big bang cosmology, in particular the origin 
of density fluctuations and the spectrum of cosmic mi- 
crowave background anisotropics. The basis of the whole 
mechanism is the stretching of quantum fluctuations gen- 
erated at sub-Hubble scales due to the exponential ex- 
pansion of the spacetime during inflation. This model, 
however, has a serious "problem" : if we consider a plain 
scalar-field-driven inflationary model — say, chaotic in- 
flation — the period of inflation lasts so long that the 
wavelengths of the fluctuations which at present corre- 
spond to cosmological scales were sub-Planckian at the 
begining of the inflationary phase. Therefore, the evolu- 
tion of fluctuations at such scales is supposed to follow 
different rules from those provided by the standard the- 
ory of cosmological perturbations — which is based on 
quantum field theory and general relativity. The set of 
rules expected to hold above this energy scale is the so- 
called trans-Planckian physics (TPP from now on). The 
TPP could lead to deviations from standard predictions 
on the cosmic microwave background radiation (CMBR), 
which probes the scales we mentioned before. The ques- 
tion to be asked is "whether the predictions of the stan- 
dard cosmology are insensitive to effects of TPP" . This 
is the precise statement of the trans-Planckian problem 

an. 

In the present paper we adopt, to mimic TPP, a series- 
expansion expression for the nonlinear dispersion relation 
uip (Eq. (TIT)|) below) , first suggested by [|[ , to modify the 
equation of motion of a test scalar field in a de Sitter 
spacetime. We argue that, in this framework, the WKB 
approximation is valid only for k 3> aH (or k p = k / 'a S> 
H with k p the physical wavenumber and H the Hubble 
factor). As a consequence, the perturbative approach 
based on the WKB approximation, used in Q to tackle 
this issue, should be used with caution. 

The outline of the paper is the following. In Sec. HTlwc 
present the standard approach for first-order cosmolog- 



ical perturbations and then we restrict ourselves to the 
case of a test scalar field in a de Sitter spacetime. We 
also introduce the WKB approximation in this section. 
In Sec. IIIII we adopt a particular way to approach TPP: 
a modification of the equation of motion due to a non- 
linear dispersion relation. The initial conditions for the 
test scalar field are set in Sec. IIV1 In the following sec- 
tions we present a numerical calculation, a 3-piecewise 
approximation (adopted in Ref. Q) and a semianalytical 
approach to solve the problem. In Sec. I Villi we inves- 
tigate a realistic model of inflation, namely power-law 
inflation, using the same nonlinear dispersion relation. 
We then conclude in Sec. IIX1 



II. STANDARD APPROACH 

The spatially flat Friedmann-Robertson- Walker 
(FRW) metric with first-order cosmological fluctuations 
is given by (ds 2 = g llv dx tl dx u ): 
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with Dij = did j 



l/3V 2 <5ij, considering the conformal 
time rj. To first-order scalar, vector and tensor perturba- 
tions evolve independently. Vector perturbations can be 
omitted because they die away kinematically. The tensor 
perturbation hij has only two physical degrees of freedom 
(polarization states) h + and /i*: 



hij = h + e 



(2) 



where and e*j are the polarization tensors having the 
following properties in the Fourier space: 

fiy = djii k eij =0, = (3) 
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4e+«=2, e?-e*«=2, e+e«-> = (4) 



which becomes for a de Sitter spacetime 



On expanding in Fourier modes we can define 



h\(t,x) 
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a(r)) (2tt) 3 / 2 
with A = +/* and where fitk will satisfy § : 



d 6 km k {ri)e* k x , (5) 



Mtfe + [k Mtfe = 
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The scalar sector, in the case of a plain scalar field 
<^>(?7, a?) = cf>(r]) + a?) driven inflationary model, can 
be reduced to the study of a single gauge-invariant scalar 
variable defined by 



a H \ 6 



(7) 



where Ti = — . This so-called Mukhanov variable @ 
obeys the following equation of motion in Fourier space: 



/4 + fe 2 - — U.* = o 



(8) 



with z = a|r. Considering an adiabatic vacuum state 
as initial condition for those perturbations, all the statis- 
tical properties are characterized by the two-point cor- 
relation function, namely by the power spectrum. The 
dimensionless power spectrum for scalar and tensor cos- 
mological fluctuations are given, respectively, by 
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(9) 



Their dependence on the mode k is defined by the spec- 
tral index in the following way: 
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dhxPq 
din k 



d\nP h 
din k 



(10) 



evaluated at a scale k <§; aH when the mode is outside the 
horizon. For n s = 1 one has a scale-invariant spectrum 
for the gauge-invariant cosmological scalar fluctuation. 

Let us now restrict ourselves to the case of a test scalar 
field in a de Sitter spacetime where a(rj) = —l/(Hrj). In 
this case, cosmological fluctuations vanish identicaly. 

If we expand the test scalar field <j>(r],x) in Fourier 
modes 



(t>{r],x) 
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d 3 k» k (r]) e lk - 



(11) 
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then the equation of motion for each mode is given by 
/4' + 2 fa)/i fc = (12 



with 
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a 



(13) 



Q 2 (r]) = k 2 



n 



2 ' 



(14) 



We can note that Eq. (fT2")) is the same equation of motion 
as that of a tensor perturbation In this simple case, 
Eq. (fl~2|) can be exactly solved. 

The two-point correlation function is given by 



(01^,5)^,^ + ^10) = 



dk sin kr 
k kr 



and the power spectrum is 
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(16) 



For supcrhorizon modes (k -C aH) the spectrum is time 
independent and scale invariant, as one can see from the 
exact solution in the limit 77 — > — . 

Equation (| 1 2|) can also be interpreted as a Schrodinger 
equation for a stationary wave function with energy 
E = u) 2 = k 2 in an effective potential V c ff(r)) = a" /a 
which is a function of the "spatial" variable 77. We might 
be tempted to solve this equation using WKB approxi- 
mation, just as it is usually done in basic quantum me- 
chanics (QM) 0. In this approximation the stationary 
solution of Eq. (I12p is given by 
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exp 



+i / (l(n')drf 



0(77') dr\ 



(17) 



as long as the WKB parameter W is much smaller than 
1: 



W 



1 

02 



2~n 



< 1 



(18) 



Figure Q] shows the behavior of W as a function of 
the conformal time r\ in a de Sitter case. According to 
this plot, the WKB approximation holds for 77 — > — 00 
(k 3> aH, subhorizon scales) but not for values close to 
zero (k <C aH, superhorizon scales), where it equals 0.125 
(i.e, not much smaller than 1) [§]. It fails exactly where 
it is supposed to: near the classical turning point (?7tp = 
—\p2jk) and where the effective potential is too steep 
(rj — > 0~). For the subhorizon scales, where E 3> V e s, the 
solution of Eqs. (|12I14|) is a plane wave in conformal time 
with comoving frequency f2 ~ k, Eq. (|17p . as expected. 



III. TRANS-PL ANCKIAN BEHAVIOR: 
NON-LINEAR DISPERSION RELATION 

A modification of the linear dispersion relation (d.r. 
from now on) was proposed by Unruh Q for describing 



FIG. 1: Total energy (full straight horizontal line), effective 
potential (dashed blue curve), and WKB parameter (full red 
curve) as a function of the conformal time r\ for a linear dis- 
persion relation. We have used k = M p i and r\ is in units of 
l/M pl . 



FIG. 2: Nonlinear dispersion relation as a function of the 
physical wavenumber k p for a case with three different solu- 
tions (k p i, k P 2 and k p s) of the turning-point equation lo p = 
2H 2 . H = 0.5M p i and k p is in units of M p i. 



high-energy Physics in the black-hole radiation emission. 
He was inspired by sound waves, for which a linear d.r. 
ceases to be valid when the wavelength gets closer to 
or smaller than the lattice spacing. Jacobson and Corley 
[To! (see also 0) also proposed nonlinear terms in the d.r. 
that could be justified by the inclusion of higher-order 
derivatives in the Lagrangian. Other changes can also 
be introduced by arguing that the spacetime symmetries 
migh not survive at high energies 

Following this approach to mimic TPP one writes the 
comoving frequency uj as lu — a(rj)LV p (kp) — where lo p 
is the physical frequency — and assumes that lo p is a 
nonlinear function of k p which differs from the standard 
(linear) one only for physical wavelengths closer to or 
smaller than the Planck scale. Note that this replacement 
would be innocuous if the dispersion relation was linear 
in k p . 

In this paper we focus on the d.r. 



u$(kp) = k*-ak* + !3kl 



(19) 



with a > and j3 > 0, proposed in Ref. @. See Fig. [2] 
for a sketch of this function. The above expression can 
be seen as a mere series expansion, but it is also found in 
solid state physics and describes the rotons [l2j]. There 
has been suggestions fl3j for including a 3-order term in 
the nonlinear d.r. above, coming from an effective-field- 
theory approach. Since such odd-power terms violate CP, 
we do not consider them here. 

Following the discussion in Section[TTl we plot in Fig. [3] 
the analogous quantities to W (fTS)) , the total energy E = 
k 2 , and the effective potential, the latter being given by 




Vcfffa) = (ak 4 )r, 2 - (/3k 6 ) V 4 



(20) 



FIG. 3: Total energy (full straight horizontal line), effective 
potential (dashed blu curve), and WKB parameter (full red 
curve) as a function of the conformal time rj for the nonlin- 
ear dispersion relation Eq. (|19p . We have used a = 0.1325, 
f3 — 0.004375 and H = 0.5M P ; so that there are 3 turning 
points (indicated by the vertical lines). The WKB curve was 
multiplied by 2 ■ 10 -2 so its behavior in the region close to 
the turning points can be clearly seen. Note, however, that 
the WKB approximation is valid only in the region 77 — > — 00 
since W(j] 0~) = 0.125. We have used k — M p i and r\ is in 
units of 1/Mpi. 



where a = aH 2 and (3 = (3H 4 are dimensionless quan- 
tities. As in the linear case, the WKB approximation is 
valid only when 77 — > — 00. 

To compare our results with Ref. [i| we consider only 
a particular range of values for a and /3, necessary for a 
positive nonlinear dispersion relation with three real dis- 
tinct solutions for the classical turning-points equation 
lu 2 = 2H 2 (see, also, }14|). The request to have three real 
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turning-points and a nonlinear u p that differs from the 
standard linear one only for physical wavelengths closer 
to or smaller than the Planck scale has as consequence 
that H should be comparable with M p i. In fact, every- 
where in the paper we consider H = 0.5M p i. 

Defining z = 3/3 /a 2 one obtains the following ranges 
for z and a: 



3/4 < z < 1 , g(z) <a< f(z) 



with 



— - 1 - (1 - zf' 2 
2 V ' 



| - 1 + (1 - zf' 2 



(21) 



(22) 



(23) 



where z > 3/4 is the physical condition to have a positive 
nonlinear d.r.. The problem of finding the power spec- 
trum from this particular nonlinear dispersion relation 
has been tackled in Ref. [4J. In that paper the authors 
introduced an approximate piecewise form of the nonlin- 
ear dispersion relation (see Sec lVIj) and apply the WKB 
approximation for k p such that ui 2 > 2H 2 . They argue 
that the amplitude of the effects of the nonlinearity is 
proportional to A = (k p i — k P 2)/k p \ where k p i and k P 2 
are, respectively, the first and the second turning points 
from the right of Fig. [5] Still according to Ref. [4(, A 
measures the "time" spent in the region where the WKB 
approximation is not satisfied, but it actually is simply a 
measure of the distance between the two largest turning 
points — k p i and k P 2 in Fig. [21 It is clearly seen in Fig. [3] 
that A is much smaller that such time interval. 

We show below the full result (i.e, the exact solution 
within the three-piecewise approximation) and two other 
different approaches — one numerical and the other semi- 
analytical — for evaluating the nonlinear effects on the 
power spectrum. 

But before proceeding to proving the claims above, we 
shall first determine the initial conditions we assume. 



IV. THE EQUATION OF MOTION: INITIAL 
CONDITIONS 



The equation of motion for the mode function of a 
test scalar field, following the approach just introduced 
to consider the TPP effects, is given by 



2 2 

a ljz 



(24) 



with uj p given by Eq. (jl9|) . 

For k ^> aH the WKB approximation is valid (see 
Fig. [3]), the term a 2 Lo 2 — uo 2 dominates and the mode 
function is given, chosing the adiabatic vacuum (see H|), 
by: 



/J-wkb(k,r)) = 



VMM) 



exp 



-i I U>(k,T)dT 



(25) 



To proceed with either a numerical or an analytical ap- 
proach to find the solution of Eq. (f2"4"| we have to find the 
initial conditions associated to our nonlinear dispersion 
relation in a de Sitter spacetime. We consider a initial 
time r]i and a fixed value of k (which should be much 
larger than a(rji)H). 

In this case, for this fixed value k 3> aH, the comoving 
frequency oj becomes nearly equal to \/]3k 3 /a 2 and the 
WKB vacuum is given by 



its derivative by 



1 



1 



y2/3V4fc3/2 Hri 
x exp 



(26) 
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V2/3!/4fc3/2 H V 2 V2 '/ 



x exp 



and our initial conditions by 

1 1 



Hwkb{rji) 

VwkbiVi) 



v / 2/3 1 /4 fc 3/2 Hm 

1 1 ./3V4fc3/2 



V2/3!/4fc3/2 Hr$ V2 



(27) 

(28) 
H m (29) 



To be more accurate the above equations are meaning- 
ful only when 

£•4 1,6 9 h.6 

tfW^tf-alL+fjIL-jL-plL (30) 
a a r\ a 

so one obtains the following constraints on (5: 



P»2H 2 ¥ , (3> >T , , /?»a^. (31) 



k 4 



k 2 



The WKB solution is an exact one at the infinite past 
and thus the choice of the adiabatic vacuum is somewhat 
"natural" . We also recall that all vacua prescriptions are 
equivalent up to zero-order when the WKB approxima- 
tion holds [15j. 

The choice of a different set of initial conditions at a 
given rj — rji, though, can always be seen as the out- 
come from a particular (trans-Planckian) evolution from 
a different (and perhaps more consensual) set of initial 
conditions defined at 77 = 770 < f]i- In other words, the 
choice of initial conditions is equivalent to the choice of 
the physics that takes place before the moment when 
they are set. Nevertheless, that choice does not replace 
the discussion on the physics that takes place after that 
moment, while energies above the Planck scale are still 
at play. 
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FIG. 4: Correction C a p on the power spectrum, with z vary- 
ing from 0.775 up to 0.950 with step 0.025, as a function of 
A from a full numerical calculation using the nonlinear d.r. 

msj. 
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FIG. 5: Three-piecewise linear approximation to the nonlinear 
d.r. as suggested in Ref. Q. The horizontal line is \/2H. 
Regions 1,2,3 are defined by the vertical lines and labeled 
from right to left. k v is in units of M p i. 



V. NUMERICAL APPROACH 

In this section we find the correction to the power spec- 
trum, for our test scalar field in a de Sitter spacetime with 
a nonlinear d.r., using a fully numerical approach with 
initial conditions given by Eqs. (|28|29|) . This has been 
done using a C code and the GSL library [16(, where 
we have set H = 0.5M p ; and considered a fixed value 
of k 3> a(r]i)H, stopping the evolution at a time rjf for 
which k <C a(rjf)H. We have verified that the result is in- 
dependent from the value of k, namely that, as expected, 
the power spectrum is still scale independent. 

We write the power spectrum as 



( §:) [i + cw 



(32) 



The function C a p represents the correction due to the 
nonlinear d.r., which obviously depends on the parame- 
ters a and p. 

In Fig. [4] we show the correction C a p for different val- 
ues of z (going from z — 0.775 to 0.950 with step 0.025) 
in function of A. We can clearly see that A is not the 
only parameter that plays a role in the calculation of the 
nonlinear effects on the power spectrum. Indeed, effec- 
tive potentials with the same A have different heights, 
depths and steepnesses, all of which influence effects on 
the power spectrum. Even for small A, the correction can 
be large depending on the values of a and (3. Therefore, 
it can hardly be considered a perturbation. 

/From Fig. we can note that for the smaller value 
of z one obtains a minimum for C a p at A > (which 
corresponds to a finite value of /3 for a fixed z) . Besides 
we have no correction (C a p = 0) for particular values 
of z and A. This means that we could still have no 
correction in spite of finite a and f3. In other words, 



this nonlinear d.r. yields no correction whatsoever to the 
power spectrum if the parameters happen to be around 
those values. 



VI. LINEAR APPROXIMATION 

Since the WKB factor (fl~5|) is not much smaller than 1 
except for k 3> aH , one cannot use the WKB expression 
(fl~7f for the solution of (f2"4|) (except in the forementioned 
region, of course). Nevertheless, there is indeed an ex- 
act solution of Eq. (|24p if we approximate the nonlinear 
dispersion relation by a straight lines, as done in Ref. (4j. 

If we write u>i(k p ) = Aik p + Bi for each of the regions 
i = 1,2,3 (see Fig. [5]), the solutions are 



+ d t W w 



iB i /9 



~fT' V4 ~ 



(33) 



where 7, = (Bi/H) 2 and Wm,w{' •>'•>') are Whittaker 
functions [13] • We pick Ai and Bi such that the posi- 
tions of the maximum and minimum of the piecewise d.r. 
coincide with those of the full nonlinear case and that 
the first and the second "turning points" k pl and k P 2 are 
the same as well (see Fig. [SJ- In this way we have the 
same A of the exact nonlinear d.r.. Of course, we pick 
A3 = 1 and £>3 = 0, corresponding to the linear d.r. in 
the long- wavelength (small k p ) limit. 

The coefficients ci, d\ are given by matching the 
asymptotic behavior of the previous equation to the ini- 
tial conditions yielded by the WKB approximation (see 
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below). The former is given by [Ei 



-iAi kr} 



(2iA 1 kr ] ) +l - 



e +lAlkri (2iA 1 kr])- 
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Cl 



r(i + 2i/) 
r(i + j/ + k ) 

r(i + 2i/) 
r(| + »/-/s) 



(34) 



where f = y9/4 — 71 and ft = —2B\jB. Comparing 
Eq. ([211) to the expected WKB solution 



IJ-wkbiv) 



1 




^^exp 


[-' f 




L Jrii 






1 / 



aj(rj)dT 



y/2Aik \r)i 



— e 



-iAi k(r]— 7]i ) 



we get 



Cl : 



(2iAikrji 



>.-iB 1 /H 



(35) 



(36) 



Note that the expression for d\ reduces to the expected 
form when A\ = 1 and B\ = 0, which correspond to the 
usual d.r.. 

The amplitude of the growing mode can be calculated 
as is usually done in QM. Since we know the exact so- 
lution in each region (Eq. ([33]) ). all we have to do is to 
match them and their derivatives at the boundaries. In 
matrix notation, at the boundary between regions 1 and 
2, we write: 



(37) 



where the superscripts (1,2) indicate the value of the 
subscript i in Eq. (|33[) . r/i 2 = — fc pl2 / (Hk) is the value of 77 
at that boundary (see Fig. [SJ and (') indicates derivative 
with respect to rj. We will write the above equation in a 
more compact form as 



) 






VdJ 


) 









M 1 -C 1 = M 2 -C 2 



(38) 



with an obvious notation. Analogously, we can write the 
matching at the boundary of regions 2 and 3 as 



'2:',) 



'(3), 



and as 



M a -C 2 = Ma-C 3 



) 


(c 2 \ 




\d 2 ) 






) 





(39) 



(40) 




FIG. 6: Correction to the spectrum using the 3-piecewise ap- 
proximation, with z varying from 0.775 up to 0.950 with step 
0.025, as a function of A. 



Such a compact notation allow us to write C 3 in terms of 
Ci as 



C 3 =M^ -M.-M^ 1 -Mi-d 



(41) 



The coefficient of the growing mode is given by d 3 — 
the second component of C 3 — since it is the coefficient 
of the divergent Whittakcr function when 77 — * 0~: 



w 



iB 3 V 1 - 4 73 



MA^ri 



2H 1 2 > j kr} 

while Wm — * at the same limit. Therefore, 

-i , 



(42) 



krj 



(43) 



as r\ 



. We can thus write the power spectrum as 
fc 3 no 2 
2^2 
H 2 k 



2tt 



-\d 3 



(44) 



in the limit 77 — > 0~ (fc <C aH). This is, as expected, 
scale invariant. 

In Fig. [5] we have plotted the correction of the spec- 
trum, due to this three-piecewise approximation, as a 
function of A. The reason for such large values is that 
the nonlinear d.r. can be qualitatively different from the 
linear one (u> ~ fc) at large \rj\ even for small A. 

As it is clear, comparing Fig. [6] to Fig. [4J the re- 
sults obtained with this three-piecewise approximation 
are pretty different from the numerical one obtained with 
the exact nonlinear d.r.. 



VII. SEMIANALYTICAL APPROACH 

In this section we use k p as the independent variable. 
Since it is proportional to rj (k p = k/a = —rjHk), there 
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is no particular advantage to choose either way, but this 
variable is more "physical" and thus one can rely on 
her/his physical intuition. For a linear d.r., Eq. (12]) 
is then written as 



d 2 ^ k 
dk 2 



1 2 . 



Using a nonlinear d.r. amounts to the substitution 



1 
IP 



1 ^p(k p ) 
H 2 k 2 



(45) 



(46) 



in the previous equation. Although there is no exact so- 
lution of Eq. (|4"5|) with the substitution (|4"6")) when one 
uses the d.r. (fl9|) . we still can get a fair analytical ap- 
proximation by writing 



bp 



1,2 



iy2 ;• J J 2 p ■ ■ > 



k- 2 



1 

H 2 



H 2 > 



a 1,2 

Tr^^p 



l 

ff 2 



a p , 1 2_j jL < jL 



(47) 



where c, d and the matching point k* are defined by 
requiring that 



Xj-^(/C*) (A/* ) S (A^*) 



(48) 



which correspond to K = (#7/3) 1/6 , d = -((3/H 2 ) 1/3 
and c = 2(/3/i7 2 ) 1 / 3 . The approximated expressions have 
the same limiting behaviors (both at k p — ► and fe p — > 
oo ) than the original one. That feature provides another 
piece of information on the dependence of the outcome 
on the details of the nonlinear d.r. for intermediate k p , 
as opposed to its small- and large-scale limits. See Fig. [7] 
for a qualitative comparison. 

The calculation itself is carried out by the usual ap- 
proach in QM, as in the previous section. The require- 
ments that the "wave function" and its first derivative 
are continuous at fc* are easily accomplished since the 
solutions of the equation of motion (145I46P with the ap- 
proximations (|47|) are known analytically. For fc > fc, we 
find: 



d 2 h(k p ) 



d^ 



-S?(fc p )/ 1 (fc p ) = 0, 



/i(fc p ) = ciexp 
+di exp 



i 

H 



^kl 



i 

H 



'2V0 kp+ 3 

■^k p + ^-kl 
p 3 p 



(49) 
Hh(k p ) + 

Ht 2 (k p ) (50) 




FIG. 7: Effective potential V c s{k p ) (thick black line) and its 
approximations (thin colored lines) given by V e ff(k p ) — E — 
£ 2 (fc p ) and Eqs. (gHJ, for a = 0.1325, = 0.004375 and 
H — 0.5M P ;, as a function of k p (in M p i units). The straight 
horizontal line is E = 1/ H 2 . The vertical line marks the value 
of fc* (see text for definition). 



where Htx{k p ) = HeunT(A,0, B,+i() and 
Ht2(k p ) = HeunT(A,0, B,—iQ are triconflucnt 
Heun functions [H| , with 



A = 



B = 



2H 2 y/<3 



2H 2 (3 2 



2/3 



1/3 



-1 + cH 2 



C = (3V2~p 1/A H) k p . 



1/3 



«y (51) 

(52) 
(53) 



We can constrain the coefficients (ci , d\) in Eq. ([50)1 by 
comparing the asymptotic behavior of the above equa- 
tion to the expected WKB one, Eq. (|2"6"|) . The asymp- 
totic limit of interest here is taken along a Stokes line 
(arg(iC) = t/2) of the triconfluent Heun function, which 
means that there are two equally dominant terms in the 
asymptotic expansion: 

HeunT(A, Q, B, iQ ~ + 6exp(-< 3 )] (54) 

as C — ► +oo, where a and b are constants. Since neither 
one is dominant over the other, both must be taken into 
account. The forementioned comparison yields 



d\b — c\a = 
— c\b + d\a = C 



(55) 



where C = i/[(3HV2P 1/4 ) 1/3 Vk}. There are only 2 equa- 
tions for 4 unknowns, which should be expected from the 
lack of a dominant behavior in the asymptotic expan- 
sion along a Stokes line, as mentioned above. One could 
choose whichever 2 of the above parameters (say, a and 
b) to be specified by comparing the final result to the 
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outcome of the numerical calculation using the approxi- 
mated effective potential. Formally, that determination 
is supposed to be done at every value of a and /?, which 
would render our semianalytical approach useless. Nev- 
ertheless, we have found that writing b and d\ in terms of 
a and c± yields a qualitatively good behavior for different 
values of z, i.e, for different pairs {a, /3} when compared 
to the numerical calculation (see Fig. [5] below). Such 
procedure yields 



o = —a 



di 



di 



-C ± y/C 2 + 4a 2 c 2 
2a 



(56) 



(57) 



The sign in the above equation was again numerically 
determined to be the lower one (see below). Following 
this line of reasoning, we used a = 0.1 and C\ = —1.8 for 
all z. 

The last step of the semianalytical procedure is the 
evolution in the second region, when k < k* : 



d 2 h(k p ) 



£!(fcp)/i(fcp) = 0, 



(58) 



whose exact solution is 

h{k p ) 



1 / 3 y/akl 

c 2 — W M \D,-, V 



1 



k„ Ww \ D >4' 




(59) 



where D = (1 - dH 2 )/ '(4Hs/a). 

The coefficient di determines the amplitude of the per- 
turbations since the function Ww(k p ) is the growing so- 
lution: 



lim 



Ww(k p ) _ V^H 1 ? 4 1 



2a 1 /8l\£) k p 



(60) 



where £ = ^| — \ I f^ a - J • The above expression allows us 
to write the spectrum as 



a 

kV~H 



8 7 ra 1 /4 r 2(£) 

H 
2^ 



\d 7 \ 2 



(1 + C afj ) , 



(61) 
(62) 
(63) 



in the superhorizon limit k p <C H. The coefficient c?2 
is determined by the forementioned procedure, requiring 
that the function and its derivative are continuous at 
k = fc* . In matrix notation, it can be written as 



dgi(k p ) 
dk v 



=fe» 



( dhijkp) \ 
\ dk p J 




I 92{K) 



h 2 {k*) 



(>2 



dgi{k p ) 



(64) 



where 
gi(k p ) = exp 

hi(k p ) 
92(k p ) 



exp 
1 



k + V^ fc 3 

H P 3 p 

/ 3 %/a/c 2 



h 2 {k p ) 



1 




ffti(*p)(65) 
ffi 2 (fc p )(66) 

(67) 
(68) 



which can be written in a more compact form, as in the 
previous section, as 

Mi ■ Ci = M 2 ■ C 2 . (69) 

As before, one can invert such equation and write 

C 2 = M 2 l ■ Mt -Ci, (70) 

which shows once more that d 2 is a linear combination 
of ci and d\ . Since we have fixed c\ from the beginning, 
we are not able to find analytically the dependence of 
d 2 on k and, therefore, to say if the spectrum is scale 
invariant, as it should be. Nevertheless, we approach the 
same problem numerically and we verify that the problem 
exhibits such invariance. 

In order to be able to measure the acuracy of the 
semianalytical approach, we have used the approxima- 
tion (|47p and numerically evolved the initial conditions 
(|28l[29|) . The evolution was split in two pieces: for k > k* 
and k < fe*. Such calculation, besides fixing the sign in 
Eq. ([57]) . also yields the "best" (i.e., more robust with 
respect to changes in a and 0) values of c\ = — 1.8 and 
a = 0.1. 

In Fig. [5] we show the function C a p for different values 
of z as a function of A, as given by the numerical evolu- 
tion using the approximated expressions for the effective 
potential and as given by the semianalytical approach 
with the same approximation. We can see that the cor- 
rection depends also on z, showing that A is not the only 
variable to work with. 

The semianalytical solution is not expected to work 
well for small /3, which means, for a fixed z, small A. 
For that, we should have taken into account the next-to- 
leading order terms in the asymptotic expansion of the 
Heun functions, Eq. (|54p . That would have, however, 
introduced new parameters that would have had to be 
fixed by comparing to the numerical solution once more. 
Since we see no advantage in having a large number of 
such parameters, we did not do so. 

We note that we reproduce a minimum Cap at a finite 
A as in the previous section. This approximation yields a 
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z=0.775 



z=0.95Q 



z=0.775 




z=0.950 



with rji = U/(l — p). The equation of motion for the 
Mukhanov variable, Eq.©, becomes (using r\ as the in- 
dependent variable) 



k 

a 



fJ-sk = 



(76) 



and the canonically normalized solution, associated with 
the adiabatic vaccum for fc 3> aH, is given by (see, for 
example, (2fj |) 



711 j 



1/2 



H^i—krj) , 



(77) 



with v = 3/2 + l/(p— 1). So, in the long- wavelength limit, 
we obtain the following scale dependent power spectrum 



FIG. 8: Correction C a /3 as a function of A when z varies from 
0.775 to 0.950 (as shown by the labels) for the semianalytical 
approach (dashed blue lines) and for the numerical evolution 
with approximated effective potential (|47|l (full red lines). 



correction C a p with a good qualitative behavior as com- 
pared to Fig. [4] but a large difference in the magnitude 
of the effect. 



VIII. POWER-LAW INFLATION 

A numerical approach can also be applied to find the 
correction to the power spectrum in a realistic model 
of inflation for the nonlinear dispersion relation under 
consideration. In this section we restrict ourselves to the 
case of power-law inflation [l9| where, in proper time, 
a(t) ~ t p , with p > 1. This expansion is generated by an 
exponential potential 



V(<l>)=Voexp 



X 

' M p i 



(71) 



with V = ^-p(3p- 1) and A = (l) 1 2 . The scale factor 
and the homogeneus solution of the scalar field are given 
by 



o(t) 



0(t) = 4 + M P /(2p) 1/2 log- 

H 

respectively, and the slow-roll parameters by 



ei 



£2 = M pl — - -, 



(72) 
(73) 

(74) 



where V$ = dV/d<p and so on. In conformal time ?y, the 
scale factor and the Hubble parameter become 



P-i i 



(75) 



P Q {k) 



1 



(2n) 2 



_1V Ir 2 (^)fc-^T (7 8 ) 



with a spectral index n s = 1 — Using this result, 

namely the dependence on k of the power spectrum, the 
range p < 60 is disfavored at 2a by the observation (see, 
for example, (2l|). 

Now we want to see if the introduction of the particular 
nonlinear d.r. given in Eq. (|19[) changes the dependence 
of the power spectrum on k. Let us set our background: 
Eq. ([76]) can be written, using Eq. (|75j) . as 



k z 



2p 2 - p 1 



(f — p) 2 rf 
For p 1 , this can be approximated as 
& + [k 2 - 2a 2 H 2 ] fx k = , 



0. 



(79) 



(80) 



so, in this limit and following the same reasoning as be- 
fore, to obtain a positive nonlinear d.r. that gives three 
real distinct classical turning-points {ui 2 = 2H 2 ) at the 
initial time (2o| we require that 



3/4 < z < 1 



9iz) 



< a < 



(81) 



with g(z) and f(z) given in Eqs. (|22I23|) . 

The new equation of motion will be given by 



k 2 fc 4 fc 6 
— -a— +(3-jr 



2p 2 -p 1 



(1 — p) 2 T} 2 



Mfc = 



(82) 

and, repeating the calculation in Sec. IIVI with the new 
background, one obtains the following initial conditions 



1 



V2(3 l / A k*/ 2 

1 



1 



(83) 

^fc3/2 (84) 



v / 2/3 1 /4fc3/2 p-lrji y/2 

We now proceed with the numerical analysis taking 
£j = (p — 1)/M p i (and thus rji = —1/M p i). The restric- 
tion to apply the adiabatic initial condition {—krji S> 1) 
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ln(k) 



ln(k) 



FIG. 9: Logarithm of the power spectrum in a power-law 
model of inflation in function of ln(fc) for a linear d.r.. k is in 
units of M p i. 



FIG. 10: Logarithm of the power spectrum in a power-law 
model of inflation using the nonlinear d.r. (|12[) in function of 
ln(fc). k is in units of M p i. 



becomes equivalent to k 3> M p i; therefore, we are indeed 
in the energy scale of TPP. Besides, with those initial 
conditions and in the limit p ^> 1, we have H(i]i) ~ M p i 
and the new d.r. becomes very different from the linear 
one only for k — 0(M p i). For the numerical analysis we 
take z = 0.80, a = /(0 -^+^ 2 °- 80) and the limiting value 
p = 60, making the analysis for k between 100M p ; and 
300-Mp/ and going from rji to a time for which k <C aH. 
We give in Fig. [5] the logarithm of the power spec- 
trum calculated for the standard linear case in function 
of ln(fc), while in Fig. [I0]the logarithm of the power spec- 
trum calculated for the nonlinear d.r. case in function of 
ln(fc) for many different points inside the aforementioned 
range of k. As one can see from this plot, the power 
spectrum for the nonlinear case still has a power-law de- 
pendence on fc, but this d.r. produces a change both in 
the normalization factor and in the spectral index when 
compared to the linear case. Fitting the data in Fig. [TU] 
we obtain a spectral index smaller than that of the lin- 
ear case: 0.857 instead of 0.966. The new spectral index 
would be almost the same as that of the linear case if we 
have taken p ~ 15. One obtains similar results starting 
from different values of p or with different values for the 
parameters a and f3. Therefore, the introduction of the 
nonlinear d.r. gives a stronger dependence of the power 
spectrum on k and, as a consequence, a possible different 
range of disfavored values for p. 

A different approach to this problem in a power-law 
model is shown in [22]. The authors use the minimum- 
uncertainty principle to fix the initial conditions and also 
find corrections to the power spectrum. We also would 
like to underline that our results are not in disagreement 
with Ref. [23| . Indeed, we consider a regime for which 
H is of the order of M p i while the results of [23| assume 
H < M p i. 

Using the same parameters as in Fig. I10[ we have plot- 
ted in Fig. [TT] the behavior of 1/(1 + W), where W is 




FIG. 11: Plot of the function 1/(1 + W) in terms of the 
wavenumber k (in units of M p i) and the conformal time f] us- 
ing the same parameters as in Fig. [10] (see text for the precise 
values). One can clearly see that the WKB approximation 
holds (W <C 1) only for early times. 

defined in Eq. (| 18|) . One can easily see that the WKB 
approximation does not hold at small absolute values of 
the conformal time, which stresses our very point 26]. 

IX. CONCLUSIONS 

In this paper we have shown that the nonlinear d.r. 
(j!9[) presents nonlinear results, i.e, which are not shown 
in the perturbative approaches used in previous papers 
in the subject. In particular, the correction C a p can 
be almost as large as 10 depending on the value of the 
parameters a, /3. Therefore, it can hardly be considered 
a perturbation. 

Our results agree, in part, with the literature on the 
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effects on the power spectrum obtained using different 
approachs to mimic TPP. In fact, also using different 
approachs, the effects are always of the order of H/M p i 
(see 0] and references cited therein) which in our case is 

0(1). 

We have compared previously used approximations 
and introduced a new one, which allowed a semianalyti- 
cal calculation. 

The fundamental point behind such a large correction 
factor lies in the break of the WKB approximation at 
early (conformal) times. If the d.r. at this moment is 
highly nonlinear, as it happens here, the corrections are 
bound to be large [24j . 

However, from Fig. |4] we can see that for particular 
value of z and A we obtain no correction to the power 
spectrum. Namely, that we could still have no correction 
in spite of finite a and f3. In other words, this nonlin- 
ear d.r. yields no correction whatsoever to the power 
spectrum if the parameters happen to be around those 
values. 

In the last section we have used a power-law model of 
inflation as the background in order to understand what 



features were particular to the de Sitter (exponential) ex- 
pansion. It is clear that one would always get a scale-free 
spectrum in a de Sitter background, but what would hap- 
pen to the spectral index n s in a different one? We have 
shown that this particular d.r. (fl"9|) yields a even smaller 
n s as compared to power-law inflation with a linear d.r., 
and therefore, it is strongly disfavored by observational 
results. On the other hand, this result suggests that an- 
other nonlinear d.r. may yield also a different correction 
(perhaps in the opposite direction) to the value of n s in 
a power-law background. 

We have not addressed the question on back-reactions 
which may play an important role, since the corrections 
to the power spectrum are not perturbative ones. Such 
implications are beyond the scope of this work and de- 
serve a separate study, which will be published elsewhere. 
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